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Abstract 

An intrinsic-state formalism for IBM-4 is presented. A basis of deformed 
bosons is introduced which allows the construction of a general trial wave 
function which has Wigner's spin-isospin SU(4) symmetry as a particular 
limit. Intrinsic-state calculations are compared with exact ones showing good 
agreement. 

PACS numbers: 21.60 -n, 21.60 Fw, 21.60 Ev 
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The Interacting Boson Model (IBM) was originally proposed to describe collective low- 
lying states in even-even nuclei. The model building blocks are monopolar (s) and quadrupo- 
lar (d) bosons. In the original formulation of the model (IBM-1) no distinction was made 
between neutrons and protons [|1|]. Later, connections with the nuclear shell model were 
investigated [0,0 and a new version was proposed in terms of neutron (s„, d v ) and proton 
(s n ,d n ) bosons, known as IBM-2 Q. The model has been widely applied to medium-mass 
and heavy nuclei, where neutrons and protons are filling different major shells. In lighter 
nuclei with N Z, however, neutrons and protons are in the same shell and a boson made 
of one neutron and one proton (known as a 5 boson) should be included. This version of 
the boson model, called IBM-3 [[|, is the simplest isospin invariant formulation of IBM. 
The three types of bosons (is, tt, and 8) form an isospin T = 1 triplet and correspond, 
microscopically, to spatially symmetric nucleon pairs with S = 0. In particular, the 8 boson 
corresponds to a spatially symmetric S = neutron-proton pair. A further extension of 
the IBM introduces the neutron-proton boson with T = or a boson, corresponding to a 
spatially symmetric nucleon pair with S = 1. This version is known as IBM-4 [[J and gives 
a proper description of even-even as well as odd-odd N ~ Z nuclei. 

The IBM-3 and IBM-4 are appropriate models for N ~ Z nuclei approaching the proton 
drip line. Such nuclei are studied intensively at the moment in particular with radioactive 
nuclear beams. Also, the IBM-4 is a reasonably simple, yet detailed model to study the 
competition between T = and T = 1 pairing, one of the hot topics in current-day nuclear 
structure physics. 

All versions of IBM are algebraic in nature and do not have a direct geometrical interpre- 
tation. Such interpretation can be achieved, however, by introducing an intrinsic state which 
provides a connection to geometric models such as that of Bohr and Mottelson 0. Intrinsic 
states have been proposed for IBM-1 for IBM-2 |TT]-0], and for IBM-3 |lJ,|T5[. Their 



primary use is to provide a geometric visualization of the model. In addition, a considerable 
reduction is achieved in the complexity of calculations, which leaves room for the inclusion 
of extra degrees of freedom. 



The purpose of this letter is to propose an intrinsic state for IBM-4. In the limit of 
strong isovector pairing it reduces to the intrinsic state for IBM-3; in general, it can be used 
for studying the competition between T = and T = 1 pairing in iV « Z nuclei. First, 
the mean-field formalism for IBM-4 is presented. This formalism is subsequently checked 
against the results of an exact calculation. 

The ensemble of bosons in the IBM-4 consists of isovector T = 1 and isoscalar T = 
bosons which have intrinsic spin S = and S = 1, respectively, to ensure spatial symmetry. 
The allowed spin-isospin combinations are thus (T,S) = (1,0) and (T,S) = (0,1). These, 
together with the orbital angular momenta i = 0, 2, give rise to 36 different bosons. The 
corresponding boson creation and annihilation operators are T^ttSg an< ^ lim,Tr,Sa where I 
is the orbital angular momentum, m is its projection, T is the isospin, r is its projection, S 
is the spin, and a is its projection. The operators ji m< TT,Sa = (—^-Y +T+S ~ m ~ T ~ a li-m,T-T,s~a 
are introduced for having appropriate tensor transformation properties. 

The construction of an intrinsic state requires two ingredients. First, it needs a basis of 
deformed bosons and secondly, it requires a trial wave function. The deformed bosons are 
defined in terms of the spherical ones through a unitary Hartree-Bose transformation, 

of _ \^ vpTrSo + Jl _ X pTrSa* t /-, x 

1 L p,Tr,Scr — 2.^ A im llm,TT,Sai !lm,TT,Sa ~ A £m S 'p,Tr,5<j; \ L ) 

im V 

and their hermitian conjugates. The deformation parameters A in these equations verify the 
following orthonormalization relations: 

E\p'TrS(j* \pTrSa _ r \pTrSa* \pTrSa _ r r /o\ 

A lm A lm -°pp'i 2^ A £m A l'm' "WW' \ l ) 

im P 

For convenience, a global label £ is used for spin, isospin, and their projections, £ = (TtSct). 
This new index plays the same role in IBM-4 as the isospin projection r does in the IBM- 
3 intrinsic-state formalism [|TjJ. The index p labels the different deformed bosons. The 
fundamental deformed boson has p = while excited ones have p — 1,2, ... , Y^i{^ + 1) — 1. 
If only s and d bosons are included, the maximum value of p is 5. Only the ground-state 
properties are considered here; so the superscript p is always zero and can be omitted 
henceforth. 
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The definition of the ground-state trial wave function follows Ref. ||15|| ; it is different 
depending on whether the system is even-even and odd-odd. For even-even nuclei with 
proton excess (the case of neutrons excess is obtained by interchanging N p and iV n ) the 
proposed trial wave function for the ground state has the form 

or, a s )) ee = tf Nn (8, a T , a fl )«4&S |0> , (3) 
where N n (N p ) is half the number of valence neutrons (protons) and 

A f (5, OCT, as) = (4=1t=i4=1t=-1 + aTfiT=lr=O fi T=lr=o) 

+ 5(4=i,r=i4=i CT =-i + as4=i CT =o4=i CT =o)- (4) 

The description of odd-odd nuclei is complicated even in the ground state, since in general its 
spin-isospin values are not known a priori. In iV = Z nuclei, which is the case of most inter- 
est, those values are known, being either (T, S) = (1, 0) or (T, S) = (0, 1). Correspondingly, 
two trial wave functions are proposed 

W,a T7 a s )} 00 ^ = At^(5,a T ,«s)4=i.=i |0> , (5) 

and 

|0(5, or, a s )) 00 _ 2 , a = At"""* (5, a T , a s )4=i. |0> • (6) 

Which of these two states is lower in energy depends on a delicate balance of the different 
terms in the Hamiltonian which in turn follow from the competition between T = and 
T = 1 pairing. 

In addition to the deformation parameters, three variational parameters, ar, oc$, and S, 
appear in the trial wave functions. The first two are related to isospin and spin symmetry 
breaking in the trial wave function. For deformation parameters independent of £ and for 
o>t = o>s = —1/2, the operator A^(<5, «r, ®s) corresponds to a bosonic pair scalar in spin 
and isospin. The parameter 5 measures the relative importance of isovector and isoscalar 
bosons in the ground state. In the limit of 5 = 0, the number of isoscalar bosons in the 



ground state is zero and the IBM-3 trial state is recovered [p!5[ . Another interesting limit is 
5 = ±1 which is obtained if the IBM-4 Hamiltonian has Wigner's SU(4) symmetry ||. In 
this case T = and T = 1 bosons are treated on equal footing. 

Given a general IBM-4 Hamiltonian, H, the ground-state equilibrium parameters are 
obtained by minimizing the expectation value of the Hamiltonian in the appropriate trial 
wave function flUJUJ^) . A general expression for this expectation value can be written as 

E(X, 8, a T} as) = e 66 h {8, oc T ,a s , 66) 
66 

+ E V L^M 5 ^T,a s ,^ 3 U): (7) 
6666 



where, 



e 6,6 — E ^imi64rn 2 6^mi^2Tn2' (^) 
l\m\l<2,m.<2, 



^6,6,6,6 ~~ ^ ^imi6,4m 2 6>4m 3 6,4"T'46 ^!mi^2m 2 ^3m 3 ^m4! (9) 

^l»7li£2™2^3™3^4"M 



and 



/i(<y, a T , as, 66) = %6 (0(Mr , as) | 0( ^ T , a5)) . ( 10 ) 



t t t t v (0(5,a T ,a s ) ^^^3% 0(5,a r ,a s )) 

/ 2 (<5, «T,a S , 6666) = 7-77 x | T77- \7 ■ I 11 ) 



The isospin matrix elements (|T0|,|TT) are calculated straightforwardly from a binomial 
expansion of the trial wave function. Furthermore, the parameters Se 1 m 1 £, 1 e 2 m 2 $,2 an d 
^imi6/2m 2 6,4w36Am 4 6 in equations (§,|) are defined as 

^ 1 mi6fem 2 6 = (hm 1 £ 1 \H\£ 2 m 2 £ 2 }, (12) 
^imi6,4m 2 6>-«3m36,4TO46 = 7 (^i m i6 ? ^2^26l ^ 1^3^36)^4^46) 



X \A + ^i^ 2 ^mim 2 ^66\/^ ^36^m 3 m 4 ^66' (^) 



where V stands for the two-body terms in the Hamiltonian H. It is worth noting that e is 
not diagonal in i. 

The energy (|7D depends explicitly on the deformation parameters A and implicitly on a?, 
as, and 5 through fi and f 2 . The deformation parameters A are obtained by minimizing 
the energy with the constraint of conserving the transformation norm, 

S[E(X, 5, a T , a s ) - £ E^Xjj = 0. (14) 
The Hartree-Bose equations for IBM-4 are obtained by deriving with respect to A|^, 

where the Hartree-Bose matrix, h^, is 



rn 2 



^3* \^2 

+ 2 £ Ve 1 m 1 z,e3m 3 Z3,e i m i z i ,e 2 m2t2 £:im:i J 4 ™ 4 f ' 2m2 f'2{S,a T , a s , ^3^2) ■ 



(16) 



There are six coupled equations of this form, which are solved for fixed values of «t, as, 
and S in a self-consistent way. This procedure yields the equilibrium deformation parameters 
A. The equilibrium values for the parameters ax, as, and 6 are obtained by an additional 
minimization. In fact, if the deformation parameters are independent of £, the parameters 
«t and as can be fixed to the value —1/2 which corresponds, as mentioned above, to a 



state with well-defined spin and isospin. As was shown for IBM-3 [|1J,|T^] this is a good 
approximation for N = Z nuclei. 

To test the present formalism, comparisons with exact calculations are carried out. Nu- 
merical calculations in the framework of IBM-4 are now possible |16| but still difficult. Also, 



only few dynamical limits have been studied. Here the following schematic Hamiltonian is 
considered 

H= a C 2 [SU TS {A)] + bC 2 [SU s (3)} + cC 2 [SU L (3)}, (17) 



where C^fG] stands for the quadratic Casimir operator of the algebra G. The first operator is 
an invariant of the algebra SUts{4) which is the boson equivalent of Wigner's supermultiplet 
algebra ||. It is worth noting here that, as mentioned in ||18|| , there are two alternative 
SUts{4) limits with the same eigenspectrum but different phases in the wave function. The 
results presented below are obtained within one them, which is associated to the election of 
the operator Y+ u (using the notation of Ref. f!8|l ). The use of the alternative limit, using Y~, 
gives identical results but changing sign to 5. The second operator in ( |T7j ) is an invariant of 
the SUs(3) algebra associated with the a (S — 1, T — 0) boson. (Its definition is analogous 
to that of SUs(3) considered in the L = IBM-4 of Ref. |TJ].) The last invariant in fliT| ) is 
an orbital deformation term associated with an SU(3) algebra which is scalar in spin and 
isospin. 

With this Hamiltonian three situations are studied. The first case corresponds to a ^ 
and c = 0. In Ref. |17| the competition of T = and T = 1 pairing was discussed using 
this Hamiltonian. This is a relevant test for the formalism presented here since it explores 
the spin-isospin degrees of freedom which makes the main difference of IBM-4 with respect 
to previous versions of IBM. In this case the mean-field and exact calculations are almost 
identical although the exact calculation is always slightly lower in energy. This can be 
appreciated in table | where exact and mean-field ground-state energies (in units of a) are 
given for N = 5 and N = 15 bosons for selected values of b/a. The value b = yields a 
Hamiltonian with the SU (4) symmetry, and degenerate lowest T = and T = 1 states. 
Negative values of the ratio b/a favor T = while positive values favor T = 1 pairing. 
The expectation value of the schematic Hamiltonian ([17]) with c = is independent of the 
deformation parameters. The minimum of the energy occurs for £ independent parameters 
and o>t = o>s = —\- This is so because it has no orbital dependence. It is worth noting that 
the variational wave functions (|3f||) contain for special values of 5 the lowest eigenfunctions 
of the SU T (3) ® SU S (3) limit @ (5 = oo) and of the SU TS {A) limit @ (5 = -1). 

The second case corresponds to a ^ and 6 = and includes the deformation term 
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C>2[SUl(3)]. The exact ground-state energy is known analytically: 

E = a A(A + 4) + c 2N(2N + 3). (18) 

This is so because, for sufficiently large negative c, the ground state belongs to the SUts{4) 
representation (0, A, 0) and the SUl(3) representation (2N, 0), where N is the boson number 
and A = T for even-even nuclei and A = 1 for odd-odd N = Z nuclei [jTj| . The corresponding 
calculation is also performed with the mean-field formalism presented here and the exact 
results are reproduced. The calculation gives an intrinsic state with (3^ = Xio/^oo = \/2, 
ax = as — ~ \ an d 5 = — 1 which is an eigenstate of the Hamiltonian ([17]) with a / and 



b = 0. This result is similar to that obtained for the intrinsic state of IBM-1 [|i9 |. 

The last case considered is the general one with a, b, c different from zero. The ground 
state still belongs to the SUl(3) representation (2N, 0) and hence the contribution to the 
ground-state energy coming from c C2[SUl(3)} is diagonal. The other two terms in the 
Hamiltonian can be diagonalized as in Ref . ]T7[] . Thus the exact energies are those calculated 
in table | (under 'exact') plus c 2iV(2iV + 3). This calculation is repeated with the mean-field 
formalism and produces an intrinsic state with the same SUl{3) symmetry as the exact one, 
(2N, 0). As in the preceding case = A20/A00 — V%, &t = «s = —\, but now 5 7^ — 1. The 
mean-field energies are those given in table | (under 'mean field') plus c 2N(2N + 3). 

These results demonstrate that the present mean-field formalism has the correct in- 
gredients for reproducing the full IBM-4 calculation. In addition, this formalism allows 
calculations for an arbitrary number of bosons and a general Hamiltonian, not necessarily 
corresponding to a dynamical symmetry limit of the model. 

To summarize, a Hartree-Bose mean-field approximation for IBM-4 has been presented, 
along with trial wave functions valid for even-even and odd-odd nuclei with N = Z. The 
trial wave functions include boson correlations in the spin and isospin spaces. Comparisons 
with exact calculations show good agreement from which can be inferred that the present 
formalism gives a good approximation to the full diagonalization. The aim is now to consider 
more realistic IBM-4 Hamiltonians that include both types of pairing (T = and T = 1), 



and a spin-orbit coupling as well as more general quadrupole deformation terms. This will 
enable the study of the interplay between single-particle, spin-isospin, and orbital degrees 
of freedom. This work is currently in progress. 

We wish to thank S. Pittel and D.D. Warner for useful comments. This work has been 
supported by the Spanish DGICYT under contract Nos. PB98-1111 and PB95-0123 and a 
DGICYT-IN2P3 agreement. 



9 



REFERENCES 

[1] F. Iachello and A. Arima, "The interacting boson model", Cambridge University Press, 
Cambridge, (1987). 

[2] T. Otsuka, A. Arima, F. Iachello, and I. Talmi, Phys. Lett. B 76, 139, (1978). 
[3] F. Iachello and I. Talmi, Rev. Mod. Phys. 59, 339, (1987). 
[4] J. P. Elliott and A.P White, Phys. Lett. B 97, 169, (1980). 
[5] J. P. Elliott and J.A. Evans, Phys. Lett. B 101, 216, (1981). 

[6] A. Bohr and B.R. Mottelson, "Nuclear Structure II. Nuclear Deformations" , Benjamin, 
New York, (1975). 

[7] J.N. Ginocchio and M.W. Kirson, Nucl. Phys. A 350, 31, (1980). 

[8] J.N. Ginocchio and M.W. Kirson, Phys. Rev. Lett. 44, 1744, (1980). 

[9] A.E.L. Dieperink and O. Scholten, Nucl. Phys. A 346, 125, (1980). 

[10] J. Dukelsky, G.G. Dussel, R.P.J. Perazzo, S.L. Reich, and H.M. Sofia, JVucZ. Phys. A 
425,93 (1984). 

[11] A. Leviatan and M.W. Kirson, Ann. Phys. 201, 13, (1990). 

[12] J.N. Ginocchio and A. Leviatan, Ann. Phys. 216, 152, (1992). 

[13] S. Pittel and J. Dukelsky, Phys. Rev. C 32, 335, (1985). 

[14] J.N. Ginocchio and A. Leviatan, Phys. Rev. Lett. 73, 1903, (1994). 

[15] J.E. Garcia-Ramos, J.M. Arias, J. Dukelsky, E. Moya de Guerra, and P. Van Isacker, 
Phys. Rev. C 57, R479, (1998). 

[16] O. Juillet, unpublished. 

[17] P. Van Isacker and D.D. Warner, Phys. Rev. Lett. 78, 3266, (1997). 

10 



[18] P. Van Isacker, J. Dukelsky, S. Pittel, and O. Juillet, J. Phys. G 24, 1261, (1998). 
[19] A. Leviatan, Phys. Rev. Lett. 77, 818, (1996). 



11 



TABLES 

TABLE I. Exact and mean-field energies of ground states, and their isospins, for selected values 
of b/a. The cases shown correspond to N = Z odd-odd nuclei with N = 5 and N = 15 bosons, 
respectively. 
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